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Abstract: The Monte Carlo generator Major 1.5 simulates the production and 
decay of heavy Majorana neutrinos via lepton mixing or exchange of 'light' right- 
handed W-bosons in deep inelastic scattering, i.e. e ± p — > NX — ► e ± W T X or v e ZX. 
Physics and programming aspects are described in this manual. 



1 Physics of included processes 

The standard model of electroweak and strong interactions has been remarkably successful in 
describing the experimental data. Nevertheless, it cannot be the ultimate theory due to its 
theoretical shortcomings, e.g. for the understanding of the mass and family structure of quarks 
and leptons and the many free parameters. Attempts to solve the theoretical problems of the 
standard model have been made based on various theoretical grounds and several extended 
theories have been suggested and studied in detail. These theories are normally based on some 
larger symmetry group which unifies the interactions and is spontaneously broken down to the 
standard model gauge group. Particular attention has been given to models with an additional 
£7(1) symmetry or with left-right symmetry in the form SU(2)l <8> SU(2)r £g> U(1)b-l, which 
can both be subgroups of the unification groups 5*0(10) and E$. The questions of neutrino 
masses and lepton number violation, which are of general importance (for a review see [|TJ), 
enter explicitly in this context and, in particular, heavy Majorana neutrinos may be present 
i- 

The simplest process, in terms of minimal new physics assumptions, to produce heavy 
Majorana neutrinos in ep collisions is through normal charged current interactions as illustrated 
in Fig. [I]. The cross-section is suppressed, not only by limited phase space due to the large 
neutrino mass m^r, but also by the required mixing £ between light and heavy Majorana 
neutrinos [§, [3|. The light Majorana neutrinos can be identified with the normal neutrinos 
u e , and u T , whose small masses is naturally explained through the see-saw mechanism. The 
mixing in the leptonic part of the charged current interaction is given by a unitary Kobayashi- 
Maskawa type matrix V and the matrix £. A power series in £ gives the connection between 
the weak eigenstates vl, vr and the Majorana mass eigenstates is, N 0, i.e. 

u L = ±=£*{y + £N +...), v R = i±pHA - + •••)■ (1) 
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Figure 1: Production (in ep collisions) and decays of heavy Majorana neutrinos. 



The heavy Majorana neutrinos will decay into either veZ°, £ + W~ or £~W + where the weak 
bosons can be taken as being on-shell since we only consider Majorana neutrino masses above 
the boson masses (mjv > tub) due to mass limits Using the narrow width approximation for 
the heavy Majorana neutrino propagator in Fig. [TJ, the cross-section for the process e~p — > NX 
where N -> £±W* or JV -> z/ £ Z° is 
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For incoming positron the last row in Eq.(^j) is replaced by 



x jC • d(x, jj 2 ) + s(x, jj 2 ) +D- u(x, jj 2 ) + c(x, jj 2 ) | 



(3) 



Another possibility to produce right-handed neutrinos is through the exchange of right- 
handed PU-bosons, Wr, see e.g. ||. If it is assumed that the Wr is lighter than the heavy 
neutrino and that there is no mixing with standard model gauge bosons then the same cross- 
section formula can be used except that G 2 F mw\{V£,) e N\ 2 = giK^OeJvl 2 * s replaced by the 
right handed gauge coupling, g R , and everywhere else m w is replaced by m WR . To simplify 
the picture one may also assume that the heavy neutrino only decays according to the charged 
decay modes, A" — > £ Wr, thus avoiding any extra assumptions on the Z'-mass. Finally, if 
only two-body decays are considered, then Wr only decays into quark-antiquark pairs due to 
the larger mass of the right-handed neutrino. It is difficult to obtain model-independent limits 
on irtw R and according to || the Wr can be as light as the ordinary Wl if extreme fine-tuning 
is allowed. 

The independent variables are chosen as the normal deep inelastic scaling variables x and y 
together with the transverse momentum p?± and rapidity yi of the final state lepton £ from the 
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Majorana neutrino decay. With four-momenta as in Fig. 1, s — (P + k) 2 , s = xs, W 2 = (P + q) 2 
and Q 2 = —q 2 . As usual, x = Q 2 /2P ■ q and y = P • q/P ■ k, but they are globally constrained 

2 2 

to < x < 1 and < y < 1 #. The phase space is given by 

- s (ptl + rn™ r lf < x <1 (4) 
s — m 2 N — J(s — m^) 2 — 4 sp™™ 2 s — m 2 N + J{s — m 2 N ) 2 — 4sp™± 2 

< y < — ^ (5) 



2s ~ ~ 2s 

.2 ™2 



o< m <^PV3 (6) 



2m^ 



m 2 N -m 2 B - 2p i± p N± m 2 N -m 2 B + 2p i± p N± 

< cosh(?/jv - ye) < = (7) 

2m N ±p e± 2m N ±pe± 



The indices iV and denote the Majorana neutrino and the W boson exchanged in the 
t-channel, whereas B and i denote the weak boson and lepton from the Majorana neutrino 
decay, respectively. The symbols ye, y^, m, m± - 



^m 2 + p\ and p± refer to rapidity, mass, 
transverse mass and transverse momentum, respectively. All frame dependent quantities are 
understood to be in the lab frame. is the total width of the Majorana neutrino and the 
functions u, c, d and s are the parton density functions in the proton evaluated at a suitable 
scale fi. 

The coefficient functions A, B depend on the decay products of the Majorana neutrino and 
are for the charged decays given by J| 
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For incoming positron the coefficient functions C, D are related to A, B in the following way, 
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(14) 
(15) 
(16) 
(17) 
(18) 
(19) 



The width of the heavy Majorana neutrino is given by 

T N = [ 2T ( N ^W*) + T(N -> u t Z) 
i 
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(20) 



Assuming the Kobayashi-Maskawa type matrix V for the lepton sector to be diagonal, so that 

BR 
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Figure 2: Branching ratios for N — ► (full curve) and N — > v^Z (dashed curve) as a 

function of the Majorana neutrino mass tojv- 

|(V£)^v| = |£i//jv|j we obtain the branching ratios shown in Fig. ||. Although all three are 1/3 
at large Majorana neutrino masses, for masses not far above the W and Z masses there is a 
substantial phase space suppression of the decay into the Z. 

For a detailed study on present limits on Majorana neutrino masses and mixings, as well as 
their phenomenology and suitable search strategies at ep colliders we refer to M. 



2 The Monte Carlo implementation 

The importance sampling method is used to generate phase space points according to the 
complete differential cross-section formula given by Eq. (§). In short, the importance sampling 
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method can be described variable transformation 



{x,y,pt±,yt) -> (H x (x),H y (y),H Pe± (p e± ),Hy e (y i )) (21) 

where H x (x) = J x h x (x')dx' , etc. and /i are the so called weighting functions. A phase-space 
point, (x,y,p£j_,yi) is first chosen from the inverse of the integrals, e.g. 

X H x [Hx,min R(Hx,rna,x -^z,min)] (22) 

where H xrain is short for H x (x min ), etc. and R is a random number, R e]0, 1[. The probability 
for keeping the point is then given by the '^-function', 

g(x,y,p e ±,yt = rrrrTu — ? — ui ? ^ 

h x {x)h y (y) h pe± (p e ± ) ^ ( 

do - 

divided by its maximum. The weighting functions /t should be chosen to mimic the characteristic 
behaviour of the differential cross-section. This makes the function g smooth which is important 
to make the rejection technique efficient. For details on the Monte Carlo methods used, see 0. 

From the phase space point (x,y,pt±,yi), the four-momenta of the particles (partons) in 
the Feynman diagram, Fig. 1, are calculated. The Lund Monte Carlo programs Lepto 6.5 
and Jetset 7.4 are then used to produce a complete final state of observable particles. The 
on-shell W/Z (from the heavy neutrino decay) is decayed with the proper branching ratios into 
a lepton pair or a qq' pair, where the polarisation of the W/Z is not taken into account. (In 
the case of a right-handed W, it is assumed to decay only hadronically since the right-handed 
neutrino is assumed to be heavier than the right-handed W.) For W/Z— > qq', parton showers 
are included to account for QCD radiation of additional partons, and this parton system is 
then hadronised using the Lund string model |T0| , |9|. Similarly, the quark coming into and 
leaving the deep inelastic scattering may radiate partons through initial and final state parton 
showers. Together with the proton remnant spectator, this parton system is hadronised with 
the Lund model. Thus, the complete 'history' of the event is generated resulting in a complete 
final state. For a more complete description of the implementation, see 0. 



3 Description of program components 

The program is written in FORTRAN 77 and consists of a set of subroutines that must be 
activated by the users main steering program, which should call the subroutine MAINIT to 
initialize the generator and then call the subroutine MAJOR for each new event to be generated. 
All subroutine and common-block names start with MA to indicate origin and avoid name 
clashes. The only exception is the real function AMGFUN which starts with AM to follow the 
FORTRAN name convention. 



3.1 Subroutines and functions 



The following subroutines should be called by the user: 
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SUBROUTINE MAINIT 

Purpose: Initiate constants and starting values. Calculate the maxima of the g-function, 

Eq. Q23D, if not given by the user. Has to be called once before MAJOR 
is called. 

Called by: User 

Calls to: LINIT, MADEFM, MADMAX 

SUBROUTINE MAJOR 

Purpose: Administer the generation of one event, calculate the four momenta and fill the 

event record. Has to be called once for every event to be generated. 

Called by: User 

Calls to: MAGENE, MAFLAV, MAERRM, LSHOWR, LUSHOW, LUPREP, LUEDIT, 

LUDBRB, LUROBO, LUEXEC 
Functions used: ULMASS, RLU, ULANGL, LUCOMP, PLU 

The following subroutines and functions are called internally: 

SUBROUTINE MADEFM 

Purpose: Define the Majorana neutrino in Jetset 7.4 code (KF=79). 



Called by: 
Calls to: 



MAINIT 



SUBROUTINE MADMAX 

Purpose: Set appropriate starting values for MINUIT |IT| and call MINUIT for the 

calculation of the maxima of the ^-function. 
Called by: MAINIT 
Calls to: MAINEW 

SUBROUTINE MASIGX(NPAR,DERIV,DIFSIG,XF,NFLAG) 



Purpose: 



Called by: 
Calls to: 



Is called by MINUIT and is used as an interface to the ^-function AMGFUN. 
Gives minus the ^-function if within range and otherwise it gives the square 
of the distance to the limit. The search for the maximum is done along the 
singularity of the ^-function (only x and y varies freely while pe± and y# are 
calculated from x and y). Stores the maxima found by MINUIT in the array 
GMAX. For details, see [0. 
MINUIT-subroutines 
MAERRM 



Functions used: AMGFUN 
SUBROUTINE MAGENE 

Purpose: Generate a phase space point according to the differential cross-section. 

Called by: MAJOR 
Calls to: 

Functions used: RLU, AMGFUN 
REAL FUNCTION AMGFUN(C) 

Purpose: Calculate the value of the ^-function, Eq. (|23|) , in a given phase space point. 

Called by: MASIGX, MAGENE 

Calls to: PYSTFU 
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SUBROUTINE MAFLAV(NINITQ,NSCATQ) 



Purpose: 

Called by: 
Calls to: 
Functions used: 



Choose flavour of initial quark and scattered quark at the boson vertex ac- 
cording to the relative cross-sections and the CKM matrix. 
MAJOR 

RLU, LUCOMP 



SUBROUTINE MAERRM(CHERRM,ERRVAL) 

Purpose: Handle errors and write error messages. 

Called by: MAJOR, MASIGX, MAGENE, AMGFUN 

Calls to: 

The following subroutines are modified versions of MINUIT subroutines from LEPTO: MACMND, 
MAIDAT, MAINEW, MAINTO, MAPRIN, MARAZZ and MASIMP. In addition the following 
subroutines (S) and functions (F) in Lepto 6.5 (L) and Jetset 7.4 (J) are used: 



Routine 
LINIT (S/L) 
LSHOWR (S/L) 
LYSTFU (S/L) 
LUDBRB (S/J) 
LUROBO (S/J) 
LUSHOW (S/J) 
LUPREP (S/J) 
LUEDIT (S/J) 
LUEXEC (S/J) 
ULMASS (F/J) 
RLU (F/J) 
ULANGL (F/J) 
LUCOMP (F/J) 
PLU (F/J) 



Purpose 

Initialize Lepto 

Include parton cascades 

Give the parton distribution functions 

Perform rotation and boost in double precision 

Perform rotation and boost in single precision 

Generate timelike parton showers 

Rearrange parton shower end products 

Exclude unstable or undetectable jets/particles from the event record 
Administrate the fragmentation and decay chain 
Give the mass of a parton/particle 

Generate a (pseudo)random number uniformly in < RLU < 1 
Calculate the angle from the x and y coordinates 
Give the compressed parton/particle code 
Provide various real-valued event data 



3.2 Common blocks 

The common-block mainly intended for communication with the program is MAUSER, which 
sets switches, parameters and cuts. In addition, all kinematic variables for a given event can be 
found in MAKINE and LUJETS in Jetset 7.4 is used to store the event record. All variables 
are given sensible default values in the block data MAD ATA and all variable names obey the 
following name convention: integers start with I-N, single precision reals start with A-C,E-H 
and O-Z and double precision reals start with D. 

COMMON /MAUSER/ MAFLAG(20),CUTM(12),PARM(30) 

Switches for controlling the program: 
MAFLAG(l) (D=9) choice of parameterisation for the parton densities in the proton. 

Is transferred to LST(15) in Lepto || which in version 6.5 gives the fol- 



lowing choices of relevance. (For information on how to use PDFLIB [12 
see LST(16) in Lepto. The parton densities are obtained from subroutine 
PYSTFU in Pythia 5.7 0.) 
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MAFLAG(2) 

MAFLAG(3) 
MAFLAG(4) 
MAFLAG(5) 
MAFLAG(6) 
MAFLAG(7) 



MAFLAG(8) 



=0: 

=1: 
=2: 
=3: 
=4: 
=5: 
=6: 
=7: 
=8: 
=9: 
= 10: 

=0: 
= 1: 
=2: 

=0: 
=1: 

=0: 
= 1: 

=0: 
= 1: 

=0: 
=5: 

=0: 
=1: 
=2: 
=3: 
=10: 
=11: 
= 12: 
=13: 



=1: 



parton density choice and parameters are controlled directly through pa- 
rameters in Pythia 5.7. 
Eichten-Hinchliffe-Lane-Quigg set 1 HI 
Eichten-Hinchliffe-Lane-Quigg set 2 



0- 



Duke-Owens set 1 
Duke-Owens set 2 



m 



MAFLAG(9) 



CTEQ2M (best MS fit) [15j. 
CTEQ2MS (singular at small-x 
CTEQ2MF (flat at small-x) p 
CTEQ2ML (large A)) 0. 
CTEQ2L (best leading order fit) |J|]. 
CTEQ2D (best DIS fit) 

(D=0) choice of factorisation scale u 2 used in the structure-functions 

Q 2 

p%j_, the squared transverse momentum of the neutrino 
the constant given in PARM(4) 

(D=0) regulates the final-state lepton status in the event record, K(I,1) 
active final-state lepton, i.e. K(I,1) = 1 
inactive final-state lepton, i.e. K(I,1)=21 

(D=l) regulates the direction of the x-axis (z-axis in proton direction) 
the x-axis along p^± 

event rotated randomly in azimuthal angle between and 27r 
(D=l) regulates the hadronisation 
hadronisation off 
hadronisation on 

(D=0) regulates the amount of written output, 
all output except MINUIT-output 
all output 

(D=0) regulates the choice of the process number (NPROC) 

I, 2 and 3 mixed according to their respective branching ratios 
e~p -> NX -> e + W-X 

e~p — > NX — > e~W + X 
e~p -> NX -> vZX 

II, 12 and 13 mixed according to their respective branching ratios 
e + p — > NX —> e~W + X 

e + p -> AX -> e+iy-X 
e+p — > AX — > i/ZX 

(D=2) simulation of QCD effects in the scattered quark and proton remnant 
system (for 2-5 cf. LST(8) in Lepto §) 

final state radiation from the scattered quark only, using a simplified treat- 
ment (LUSHOW is used on the quark and diquark system). Only valence 
quarks, z.e.u-quarks, are considered in the proton, 
initial and final state radiation 
initial state radiation 
final state radiation 

QCD switched off but target remnant as 2, 3 and 4 
error flag 
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MAFLAG(IO) 



MAFLAG(ll) 



MAFLAG(12) 
MAFLAG(13) 



MAFLAG(14) 
MAFLAG(15) 



MAFLAG(16) 
MAFLAG(17) 
MAFLAG(18) 

MAFLAG(19) 



=1: maximum of g-function violated. This is serious if the violation is large 
and in that case the user has to either set the maxima by hand (see 
MAFLAG(19)) or change the /i-functions to make the (^-function smoother 
(D=l) Error handling 

=0: no warnings and execution not stopped on error 

=1: warnings printed but execution not stopped 

=2: warnings printed and execution stopped 

(D=2) choice of /i x -function (see Section 2 for definition and Table p] for 
suitable choices) 



-2 
-1 
=0 
= 1 
--2 

--{): 
= 1: 



--2 
--1 
=0 
= 1 
=2 

=0: 
= 1: 



h x ( 
K( 
h x ( 
h x { 
Ki 
(D 



1/x 2 
1/x 
const. 
exjp(-A x x) 
xexp(—A x x 2 ) 
=0) regulates calculation of A 2 
is calculated by the program 



A x is given in PARM(8) 

(D=l) choice of /i y -function (see Section 2 for definition and Table [i] for 
suitable choices) 



h v(y) 

h y {y) 

h y (y) 

h y (y) 
h y (y) 



W 
i/y 

const. 

exp(-A y y) 

yexp(-A y y 2 ) 



--2 
--1 
=0 
= 1 
=2 

=0: 
= 1: 

=0: 
= 1: 



(D=0) regulates calculation of A y 
A y is calculated by the program 
A y is given in PARM(9) 

(D=2) choice of h Pe± -function (see Section 2 for definition and Table |l] for 
suitable choices) 
hpt X (pe±) = l/pe± 2 
h Pe± (pe±) = 1/Pt± 
h Pe± {pe±) = const. 
h n „: (pe±) = exp(-A Pt± p ei _) 

(Pe±) = Pi± exp(-A 



"PiX 

gives h 



PC. 



Pil. 2 ) 



(D=0) regulates calculation of A Pe± for — > e W T 

A pe± is calculated by the program 

A Pe± is given in PARM(10) 

(D=0) regulates calculation of A Pe± for — > uZ 

A pe± is calculated by the program 

A r , 



,„ x is given in PARM(ll) 
(D=l) choice between fixed and varying (according to a Breit-Wigner dis- 
tribution) boson-masses in the neutrino decay 
=0: fixed boson-masses 
= 1: varying boson- masses 

(D=0) regulates calculation of ^-function (see section 2) maxima 
=0: maxima of (^-functions are calculated 
=1: maxima of (^-functions given in PARM(12) to PARM(17) 
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MAFLAG(20) 



=0: 
=1: 



CUTM(1 
CUTM(2 
CUTM(3 
CUTM(4 
CUTM(5 
CUTM(6 
CUTM(7 
CUTM(8 
CUTM(9 
CUTM(IO) 
CUTM(ll) 
CUTM(12) 

PARM(l) 
PARM(2) 
PARM(3) 
PARM(4) 
PARM(5) 



PARM(6): 



PARM(7) 

PARM(8) 

PARM(9) 

PARM(IO) 

PARM(ll) 

PARM(12) 

PARM(13) 

PARM(14) 

PARM(15) 

PARM(16) 

PARM(17) 

PARM(18) 

PARM(19) 

PARM(20) 

PARM(21) 

PARM(22) 

PARM(23) 



1(^0. 



AT 



(D=0) choice between left and right handed W exchange 
left handed W exchange 
right handed W exchange 

Cuts on kinematic variables defined in Section 1: 
^min (D=0.) 

■^max (D 1.) 
2/min (D=0.) 
2/max (D = l.) 

Q 2 mm (D=4. GeV 2 ) 
QLx (D=10 8 GeV 2 ) 
Wl in (D=9. GeV 2 ) 
< ax (D=10 8 GeV 2 ) 
ftlmin (D=0. GeV) 
mmax (D=1000 GeV) 

Vimin (D=-100) 
Z/£max (D = 100) 

Parameters for input (1-20) and output (21-30): 
(D=100 GeV) mass of heavy Majorana neutrino 
(D=820 GeV) proton-momentum 
(D=30 GeV) electron-momentum 

(D=1000 GeV 2 ) fixed factorisation scale a 2 used if MAFLAG(2)=2 
(D=0.01) degree of mixing between light and heavy Majorana neutrinos 

tween 
assumed to be 1. 

(D=0.001) cut used around the singularity of the g-function which divides the 

^-function into two parts (see ALSCUT in for details). A larger cut will speed 

up the simulation but at the same time increase the risk of violating the maximum 

of the g-function and cut away a larger part of the phase-space. 

(D=l.l) safety factor, multiplies the maxima of the (yf-function. 

A x (see also MAFLAG(12)) 

A y (see also MAFLAG(14)) 

A pt± for N -> e±W^ (see also MAFLAG(16)) 

A nx for N uZ (see also MAFLAG(17)) 

Max of first part of ^-function for process nr 1 or 11 (see also MAFLAG(19)) 
Max of second part of ^-function for process nr 1 or 11 (see also MAFLAG(19)) 
Max of first part of ^-function for process nr 2 or 12 (see also MAFLAG(19)) 
Max of second part of ^-function for process nr 2 or 12 (see also MAFLAG(19)) 
Max of first part of ^-function for process nr 3 or 13 (see also MAFLAG(19)) 
Max of second part of ^-function for process nr 3 or 13 (see also MAFLAG(19)) 
Mass of right-handed 1^-boson in GeV 
Width of right-handed IV-boson in GeV 
not used 

Estimate 1 of the cross-section in pb (including cuts) for process nr 1 or 11. 
Estimate 1 of the standard deviation in PARM(21) 

Estimate 1 of the cross-section in pb (including cuts) for process nr 2 or 12. 



2 . In the case of right handed W exchange it is interpreted as the ratio be- 
the right and left handed gauge couplings, [gnj gi) 2 which normally is 
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PARM(24) 
PARM(25) 
PARM(26) 
PARM(27) 
PARM(28) 
PARM(29) 
PARM(30) 



Estimate 1 of the standard deviation in PARM(23) 

Estimate 1 of the cross-section in pb (including cuts) for process nr 3 or 13. 

Estimate 1 of the standard deviation in PARM(25) 

Efficiency 2 in the generation for process nr 1 or 11 

Efficiency 2 in the generation for process nr 2 or 12 

Efficiency^ in the generation for process nr 3 or 13 

not used 



COMMON /MARINE/ EP,EE,S,SHAT,Q2,W2,X,Y,PNT,PNL,YN,EN,PLT,YL 

Kinematical variables in a given event of use for the user: 

EP: Proton energy in GeV 

EE: Electron energy in GeV 

S: CMS-energy squared (Mandelstam s) in GeV 2 

SHAT: s = xs 

Q2: Momentum transfer squared, Q 2 = —q 2 in GeV 2 

W2: Hadronic CMS-energy squared, W 2 = (p p + q) 2 in GeV 2 

X: Bjorken-x, x = Q 2 /2p p ■ q 

Y: Standard y- variable, y = p p ■ q/p p ■ p e 

PNT: Transverse momentum of heavy p^i. Majorana neutrino in GeV 

PNL: Longitudinal momentum of heavy Majorana neutrino in GeV 

YN: Rapidity yjv of heavy Majorana neutrino 

EN: Energy of heavy Majorana neutrino in GeV 

PLT: Transverse momentum pgj_ of final-state lepton from N decay in GeV 

YL: Rapidity ye of final-state lepton 

COMMON /LUJETS/ N,K(4000,5),P(4000,5),V(4000,5) 

The variables in the common-block LUJETS are described in the Jetset 7.4 manual ||. The 
first seven entries in the event record are as follows (in the lab frame with the z-axis in the 
proton direction): 

1. Incoming electron 

2. Incoming proton 

3. Exchanged W^-boson 

4. Heavy Majorana neutrino 

5. Incoming parton before initial shower 

6. Incoming quark at boson vertex 

7. Scattered quark at boson vertex before final shower 

The decay products from the A-decay are in line K(4,4) (final state lepton) and K(4,5) (on 
shell boson) respectively. 

Common-blocks for internal use: MACROS, MAMAMI, MAMASS, MACONS, MAGSPE. 
Lepto 6.5 common-blocks used : LEPTOU, LBOOST, LFLMIX, LPFLAG, LYPARA. 
Jetset 7.4 common-blocks used : LUDAT1, LUDAT2, LUDAT3, LUDAT4. 



x The estimate is updated for each event so it should not be used until all events have been generated. 
2 Thc efficiency is given by the number of accepted events divided by the number of tries. 
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Table 1: Examples of suitable choices of h- functions at different cms energies and for different 
heavy Majorana neutrino masses. 





HERA 


LEP®LHC 


[GeV] 


300 


1200 


m N [GeV] 


100 


300 700 


MAFLAG(ll) 


2 


-1 -2 


MAFLAG(13) 


1 


-1 -1 


MAFLAG(15) 


2 


2 2 



3.3 Update history 



Updates from version 1.1 [16| to 1.3 [[T? 



inclusion of the neutral current decay N — > v^Z 

extension and redefinition of the common-block MAUSER and some internal common- 
blocks 

inclusion of a subroutine for error handling 

correction of a small error in the differential cross-section 



from version 1.3 Jl7j to 1.5: 



inclusion of possibility to have incoming positron (see MAFLAG(7)) 

inclusion of possibility to have exchange of right-handed W's which are lighter than the 
Heavy Neutrino (see MAFLAG(20), PARM(5), PARM(18) and PARM(19)) 

update to run with Lepto 6.5, Jetset 7.4 and Pythia 5.7 (see MAFLAG(l)) 



inclusion of MINUIT 111 subroutines from Lepto 



4 Usage and availability 



Major 1.5 should be loaded together with Lepto 6.5 [g], Jetset 7.4 and Pythia 5.7 ||. The 
program is a slave system, which the user must call from his own steering program. Information 
about the program, the source code and a demonstration job can be found via the WWW on 
the MAJOR homepage |http : / /www3 . tsl . uu . se/ thep/maj or/| or from one of the authors via 
email. It is not recommended to split the program and making a library of it since this can 
give problems with the use of Lepto routines. Examples of suitable choices of /i-functions are 
given in Table |I|. 
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